Recoil correction to the bound-electron g factor in H-like atoms to all orders in olL 



V. M. Shabaev and V. A. Yerokhin 

Department of Physics, St. Petersburg State University, Oulianovskaya 1, Petrodvorets, St. Petersburg 198504, Russia 

The nuclear recoil correction to the bound-electron g factor in H-like atoms is calculated to first 
order in m/M and to all orders in aZ. The calculation is performed in the range Z = 1 — 100. A 
large contribution of terms of order iaZY" and higher is found. Even for hydrogen, the higher-order 
correction exceeds the {aZ)'^ term, while for uranium it is above the leading {ctZY' correction. 



Recent progress in high-precision measurements of the bound-electron g factor for H-like carbon ||,|^ and the 
related theoretical investigations |pHlli provide a new independent determination of the electron mass. The accuracy 
of this determination presented in |11|] is three times better than that of the accepted value for the electron mass [|2| . 
This result can be improved if the theoretical and experimental uncertainties for the g factor are reduced. From the 
experimental side, an improvement of the accuracy by an order of magnitude is anticipated in the near future, as well 
as an extension of the measurements to higher-Z systems [|l3|. Investigations of the bound-electron g factor in high-Z 
systems are of particular importance since they can provide a new determination of the fine structure constant [g],^ , 
nuclear magnetic moments and nuclear charge radii. They would also create a good possibility for testing the 
magnetic sector of QED in a strong Coulomb field. At present, the theoretical uncertainty of the bound-electron g 
factor in H-like ions is mainly determined by four factors: a numerical error in evaluation of the QED correction of 
first order in a, an error caused by employing the aZ expansion for the QED correction of second order in a, an error 
resulting from the aZ expansion of the nuclear recoil correction, and, for very heavy ions, an error due to the finite 
nuclear size correction. The main goal of this Letter is to evaluate the nuclear recoil correction to the Is g factor to all 
orders in aZ and, therefore, to eliminate one of the main sources of the uncertainty for the corresponding theoretical 
predictions. 

As is known in the nonrelativistic limit the recoil correction to the \s g factor vanishes. The leading relativistic 
recoil correction is of order {aZYm/M and was evaluated in p^ , p^ (see also and references therein). General 
formulas for the nuclear recoil effect valid to all orders in aZ were derived recently in 1^. These results were confirmed 
by Yelkhovsky ||ic|] by employing a different method. In addition, Yelkhovsky presented some arguments for the 
assertion that the recoil correction up to order {aZYm/M is completely defined by the so-called lower-order term 
(which was evaluated analytically in |y and re-derived in jl0|). As a result, in the total theoretical uncertainty for 
the g factor in C'^''" was reduced to the level of 1.2 x 10^^. This leads to improving the precision of the electron-mass 
determination by factor of two. In the present Letter, we numerically evaluate the higher-order contribution to the 
recoil correction to all orders in aZ. Our results confirm the statement of [|l^ that the expansion of this term starts 
with {aZ)^ . However, we find that the {aZ)^ behavior of the higher-order term is a result of a cancellation of terms 
of order {aZy^ and {aZ)'^ (see a discussion below), while the argumentation of px[ | does not contain any indication 
of the appearance of such terms. We also observe that for all H-like atoms in the range Z = 1 — 100 the higher-order 
term exceeds the {aZ)'^m/M contribution. In particular, for the case of carbon this term is about five times larger 
than the {aZ)'^m/M term and by factor of ten exceeds its estimation given in pO[ . 

We consider a H-like atom with a spinless nucleus that is put into the classical homogeneous magnetic field, 
Aci(r) — [H X r]/2. For simplicity, we assume that H is directed along the z axis. The energy shift of a state a to 
first order in Ti. and to first order in m/M is conveniently written as the sum of the lower-order and the higher-order 
term AE = AEi^ + AEu, where (?i = c = 1 , e < 0) 
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Here, p*^ = —iV'^, n = 
D''{u;) = -47raZa'D"=(w) 



r/r, V{r) — —aZ/r is the Coulomb potential of the nucleus, SV{x.) 



47r L r cj^r J 
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-ea • Aci(x), 
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is the transverse part of the photon propagator in the Coulomb gauge, G{uj) ~ ^„ — £„(1 — iO)] ^ is the Dirac- 

Coulomb Green function, fc^ = (a|(5F|a}, \5a) = X^n"'^^" W){'^\^yW){^a — £n)^^ , and a is a vector incorporating the 
Dirac matrices. In equation (^, the summation over the repeated indices {k = 1,2,3), which enumerate components 
of three-dimensional vectors, is implicit. The recoil correction to the bound-electron g factor is defined as = 
AE/{fioTi.mj), where fxo = |e|/(2m) is the Bohr magneton and rrij is the angular momentum projection of the state 
under consideration. 

For the Is state, the analytical evaluation of the lower-order term yields P,pO| 

(aZ)4 



(4) 



The first term in the right-hand side of this equation reproduces the result of [ |l5[|l6| ] , while the second term contributes 
to order (aZ)^ and higher. 

The higher-order term, defined by equation (|^), is represented by the sum of the Coulomb, the one-transverse- 
photon, and the two-transverse-photon contribution, AEh = AE^°^^ + AE^^ + AE^^. For the Is state, we transform 
them to the form appropriate for the numerical evaluation. 
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where 8(7/) = aZaexp (— i/r)/r + iaZ[H,(j}{yr)Ta\, H — a • p -|- (3m + V{r) is the Dirac-Coulomb Hamiltonian, 
Aij =ei~ej, and (/)(yr) = [exp (— yr)(l + yr) — l]/(yr)^. The wave-function correction \da) can be found analytically 
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by the method of generahzed virial relations for the Dirac equation The exphcit form for the component of 

I (5a), that has the same angular quantum numbers as the reference state \a), is presented in H] (only this component 
contributes to the effect under consideration). 

The numerical evaluation of the expressions (B|)-(0) was carried out similarly to our previous calculations of the 
nuclear recoil correction to the Lamb shift p9|-pi|. After integration over angles, the finite basis set method was 
used to evaluate infinite summations over the electron spectrum. Basis functions were constructed from B-splines by 
employing the procedure proposed in p2| . The integration over y was carried out numerically for equation (|^) and 
both numerically and analytically for equation (|^). We mention large numerical cancellations arising in for very 
small Z if the y integration is performed analytically. In this case, numerical integration turns out to be preferable. 
The number of B-splines used in actual calculations was varied from 70 to 110. The estimated uncertainty corresponds 
to the dependence of the results on grid parameters and the number of splines and integration points. 

The correction to the Is g factor induced by the higher-order term AEn is expressed in terms of the function 
PiaZ), 

Ag^^^iaZfPiaZ). (8) 

The corresponding numerical results are presented in Table I. It is noteworthy that the one-transverse-photon and the 
two-transverse-photon contribution separately are of the order (aZ)'^ for small Z, while their sum exhibits the (aZ)^ 
behavior. This fact is clearly demonstrated in Fig. 1, where the numerical results for the ratio AgB_/[{iTi/M){aZ)'^] 
are plotted. We also note that the one-transverse-photon contribution contains terms of order (aZ)'^ which are 
cancelled when added together. Namely, the part corresponding to the perturbation of the reference state a and the 
part corresponding to the perturbation of the electron propagator exhibit the (aZ)^ behavior, when taken separately. 
We note that, in contrast to our results, the argumentation of ||lo[] , where the same gauge is considered, does not 
indicate the appearance of terms of order lower than (aZ)^. For this reason, we can not consider the argumentation 
of Yelkhovsky, in the form it is given in ||l0|| , as complete. Fitting our numerical results for small Z to the form 
P{aZ) = C51 log(aZ) + C50 + aZ{- ■ ■) yields C51 = -5.3 ± 0.5 and C50 = -6.5 ± 1.0. The uncertainties of the 
coefficients were estimated analyzing the dependence of the results on the number of parameters in the fit and the 
number of fitting points. 

In Table 2 we present the ratios Agn/Ago and (A(/l — Ago)/ Ago, where Ago = {aZ)'^m/M is the lowest-order 
correction derived in [|l^,^. As can be seen from the table, the A^h term exceeds (A^l — Ago) for all Z in the range 
Z = 1 — 100. In the case of carbon, A^h amounts to 7.7 x 10^^^, which is ten times larger than the uncertainty ascribed 
to this term in [ p^ . However, since this correction is about ten times smaller than the current theoretical uncertainty 
due to the binding QED correction, it does not affect the electron-mass prediction of ||ll|. The higher-order recoil 
correction is more important for higher-Z systems, since it grows very rapidly when Z increases. In particular, for 
uranium the higher-order recoil correction is even above the (aZYm/M term. 

In Table 3 we present the individual contributions to the Is g factor for some H-like ions in the range Z = 6 — 92. 
An error ascribed to the Dirac point-nucleus value results from the current uncertainty of the fine structure constant, 
1/a = 137.03599976(50) |l^. The uncertainty of the finite nuclear size correction was estimated as the difference 
between the result obtained with the Fermi model of the nuclear charge distribution and with the homogeneously- 
charged-sphere model. The nuclear charge radii were taken from [|3|j24). The one-loop QED correction was taken 



from |J,|5| , where it was evaluated numerically to all orders in aZ . The a QED correction includes the existing aZ 
expansion terms for the QED correction of second order in a and the known free-QED terms of higher orders 
in a (see, e.g., [^). Its relative uncertainty was estimated as the ratio of the part of the one-loop QED correction, 
that is beyond the {aZ^ approximation, to the part, that is within the {aZY approximation. The recoil correction 
incorporates the total recoil correction of first order in m/M, calculated in this work, and the known corrections of 
orders {m/MY and a{m/M) [|l^. From the table, we conclude that for low Z the theoretical uncertainty is mainly 
determined by the numerical error of the one-loop QED correction while for high Z it results from the aZ 

expansion of the QED correction and from the finite nuclear size correction. Calculations of the QED corrections 
up to the desirable accuracy seem to be feasible in the near future if we consider recent progress in calculations of the 
corresponding corrections to the Lamb shift in H-like ions [ p5|p^ . As to the uncertainty due to the finite nuclear size 
effect, one may expect that, like in the case of the hyperfine splitting |p7| , it can be significantly reduced in a specific 
difference of the bound-electron g factor in H- and Li-like ions. 

We thank T. Beier, S. Karshenboim, J. Kluge, W. Quint, and A. Yelkhovsky for valuable discussions and N. 
Lentsman for improving the language. This work was supported in part by RFBR (Grant No. 01-02-17248), by the 
program "Russian Universities - Basic Research" (project No. 3930), and by GSI. 
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TABLE L The higher-order recoil correction to the Is g factor, expressed in terms of the function P{aZ) defined by equation 
(8). 



z 


-Pcoul 


Ptii 


Ptr2 


P 


1 


-1.11414 


100.701(2) 


-80.8200(3) 


18.769(2) 


2 


-1.09754 


53.5278(6) 


-36.98689 


15.4434(6) 


3 


-1.08183 


37.44949(6) 


-22.80837 


13.55928(6) 


4 


-1.06693 


29.24592(4) 


-15.91960 


12.25940(4) 


5 


-1.05277 


24.23027(3) 


-11.90049 


11.27702(3) 


6 


-1.03931 


20.82713(1) 


-9.29387 


10.49395(1) 


8 


-1.01429 


16.47349 


-6.15902 


9.30018(1) 


10 


-0.99161 


13.78331 


-4.37646 


8.41524(1) 


20 


-0.90647 


8.09979 


-1.22907 


5.96425(1) 


30 


-0.85834 


6.08456 


-0.41612 


4.81010(1) 


40 


-0.84048(1) 


5.09672(1) 


-0.08937 


4.16687(1) 


50 


-0.85209(1) 


4.58136(3) 


0.07843 


3.80770(3) 


60 


-0.89803(3) 


4.36161(3) 


0.18668(1) 


3.65025(4) 


70 


-0.99173(9) 


4.39148(9) 


0.27839(1) 


3.6781(1) 


80 


-1.1647(1) 


4.7153(3) 


0.38378(3) 


3.9344(3) 


90 


-1.4962(9) 


5.525(3) 


0.5459(2) 


4.575(3) 


100 


-2.228(9) 


7.48(3) 


0.883(2) 


6.14(3) 
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TABLE II. The recoil corrections to the Is g factor, expressed in terms of the leading correction, Ago = {aZ)'^m/M. The 
difference AgL — Ago is defined by the second term in equation (4) and corresponds to the deviation of the lower-order term 
from [aZYm/M. Agn is the higher-order term. 



z 


(Aol " Ago)/ Ago 


AgYij Ago 


1 


-4.437732x10"** 


7.2935(8) x 10"*^ 


2 


-1.775234x10"^ 


4.8010(2) X 10"^ 


6 


-1.599074x10"* 


8.80823(1) X 10"" 


10 


-4.449468x10"* 


3.27011(1) x 10"^ 


20 


-1.794206x10"^ 


1.85414(1) X 10"^ 


30 


-4.092524x10"^ 


5.04677(1) x 10"^ 


50 


-1.190031x10"^ 


0.184956(2) 


70 


-2.514919x10"^ 


0.49025(2) 


90 


-4.672903x10"^ 


1.296(1) 


100 


-6.261303x10"^ 


2.39(1) 



TABLE III. The individual contributions to the Is bound-electron g factor in H-like ions. 





12(-,5+ 


16o7+ 


32gl5+ 




40^^19+ 


Dirac value (point) 
Fin. nucl. size 
QED, order (a/ir) 
QED, order (a/jr)^ 
Recoil 
Total 


1.998 721 354 4 
0.000 000 000 4 
0.002 323 663 7(9) 
-0.000 003 516 2(2) 
0.000 000 087 6 
2.001 041 589 9(9) 


1.997 726 003 1 
0.000 000 001 5 
0.002 324 416(1) 
-0.000 003 517 1(4) 
0.000 000 117 
2.000 047 021(1) 


1.990 880 058 3(1) 
0.000 000 038 9 
0.002 330 920(3) 
-0.000 003 523(4) 
0.000 000 236 
1.993 208 254(5) 


1.988 447 661 3(1) 
0.000 000 070 3 
0.002 333 636(4) 
-0.000 003 525(6) 
0.000 000 239 8 
1.990 778 082(8) 


1.985 723 203 8(1) 
0.000 000 113 1(1) 
0.002 336 92(1) 
-0.000 003 528(9) 
0.000 000 297 1 
1.988 057 01(2) 






7^Ge=*'+ 








Dirac value (point) 
Fin. nucl. size 
QED, order (o/tt) 
QED, order (a/w)^ 
Recoil 
Total 


1.979 392 224 9(2) 
0.000 000 272 6(2) 
0.002 345 02(1) 
-0.000 003 533(16) 
0.000 000 332 4 
1.981 734 32(2) 


1.963 137 509 5(3) 
0.000 001 231 2(10) 
0.002 369 20(1) 
-0.000 003 55(4) 
0.000 000 426 5 
1.965 504 82(4) 


1.892 114 650(1) 
0.000 023 49(3) 
0.002 505 26(1) 
-0.000 003 61(19) 
0.000 000 783 5 
1.894 640 57(19) 


1.734 947 026(2) 
0.000 453 3(9) 
0.002 884 38(3) 
-0.000 003 7(6) 
0.000 001 723 
1.738 282 7(11) 


1.654 846 173(3) 
0.001 275 0(25) 
0.003 088 93(3) 
-0.000 003 8(9) 
0.000 002 491 
1.659 208 9(27) 
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FIG. 1. The Coulomb, the one-transverse-photon, and the two-transverse-photon contribution to the ratio AgH/[{aZ)'^m/M]. 
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